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which these experiments could detect deviations from the SM Higgs predictions. 
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I. INTRODUCTION 



The Standard Model (SM) of particle physics has provided a remarkably successful de- 
scription of electroweak data up to present energies. While the minimal SM implementation 
of electroweak symmetry breaking (EWSB) relies on a single SU(2) l doublet Higgs field, the 
dynamics of the EWSB sector have not yet been directly probed and extensions of the SM 
allow for a wide variety of extended Higgs sectors consistent with all existing data.^ Indeed, 
models that address the hierarchy problem - the extreme instability of the SM Higgs mass 
parameter to radiative corrections - contain additional fundamental or composite scalar 
particles at or below the TeV scale. 

If a Higgs-like state is discovered, the next priority will be to test the Higgs mechanism 
of mass generation by measuring its couplings to SM particles ^. Experimental sensitivity to 
the Higgs couplings comes from measurements of the Higgs production cross sections, decay 
branching fractions, and total width. Prospects for the extraction of Higgs couplings from 
experimental data have been studied for the CERN Large Hadron Collider (LHC) 
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International Linear e^e~ Collider (ILC) (for a review and references see Ref. 
photon collider [l^, [l^ , and muon collider fl5|: for a recent review of Higgs boson production 
and decay at these machines see also Ref. [16|]. These measurements can be used to test 
the consistency of the measured Higgs couplings with SM predictions, and, if a discrepancy 
is found, to constrain the possible nature of the extended model. Different models give 
rise to different patterns of Higgs coupling deviations; in other words, they occupy different 
"footprints" in the space of measurable Higgs couplings. By identifying the footprint of each 
different extended Higgs sector in the space of Higgs couplings, we can determine whether 
the models can be distinguished in principle and establish a framework to map back any 
observed pattern of Higgs coupling deviations onto the appropriate underlying model. 

Our aim is to make a complete catalog of extended Higgs sectors based on the patterns 
of coupling shifts of a single neutral Higgs state h. In this paper we limit ourselves to 
models that contain only Higgs doublets and/or singlets. This allows us to avoid tree- level 
violation of custodial SU(2) symmetry and the resulting stringent constraints from the p 



^ Models with extra Higgs doublets and/ or singlets can yield better agreement with electroweak precision 

measurements than the SM; cf. Refs [l|, 
^ Determining the nature of a new observed spin-less state that is not a Higgs is also challenging 
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parameter that arise, e.g., in models containing Higgs triplets |17|. W e also limit our study 



19| for natural flavor 



to models that obey the Glashow-Weinberg-Paschos condition 
conservation, which requires that all fermions of the same electric charge get their mass 
from exactly one Higgs doublet. This allows us to avoid tree- level flavor- changing neutral 
Higgs interactions and the resulting stringent constraints from low-energy flavor physics. 



20|, and assume that the state h is pure 



We also neglect the possibility of CP violation 
CP-even. Of course, many viable extended Higgs models exist outside these categories; 
their inclusion into the framework of Higgs coupling footprints presented here would make 
an obvious future extension of this work. 

Within the above constraints we enumerate the complete set of models that can arise 
and present formulas for the shifts in the couplings of h to SM particles relative to their 
SM values. We consider only what can be learned from one Higgs state, h; of course, 
observation of additional Higgs states (CP-even, CP-odd, or charged) or other new particles 
will complement this information. We focus on the couplings that arise from dimension-four 
operators, in particular the couplings of h to W or Z boson pairs and to fermion pairs. 
Note that the hWW and hZZ couplings are modified from their SM values by a common 
multiplicative factor in any model containing only Higgs doublets and singlets. Similarly, 
our assumption of natural flavor conservation implies that the huiUi couplings for the three 
generations of up-type quarks are modified by a common multiplicative factor; the same 
holds for down-type quarks and for charged leptons. We do not give explicit results for 
the loop-induced Higgs couplings to gluon or photon pairs or 7Z because new physics at 
higher scales can generate additional effective couplings to these final states comparable in 



strength to those induced by SM loops 2l(|^. These loop-induced couplings are nevertheless 
important because they provide experimental access to the relative signs of the dimension- 
four couplings. Similarly, higher scale physics can generate dimension-six hWW and hZZ 
operators 23, 24, 25,|26|, ^'^ ^^^^ include these effects in our analysis. 

This paper is organized as follows. In the next section we briefly review the Higgs 
couplings in the SM and introduce the notation and general framework that will be used to 
describe the extended models. We then proceed to the discussion of the extended models in 



^ Indeed, efforts to distinguish new physics running in loops involved in production and decay of a lone 
Higgs state have been made for selected models |22| . 
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Sees. IIIH llVl andlVl For eaeh model we present the couplings of h in terms of its composition 
and the vacuum expectation values (vevs) of the doublets in the model; where possible we 
also invert these relations to find explicit expressions for the model parameters in terms of 
the Higgs couplings. We identify the coupling patterns that allow the different models to 
be distinguished and specify the sets of models that cannot be distinguished based on the 
couplings of only one state. We also give expansions for the couplings near the decoupling 



limit 
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29( 1 in which the deviations of the couplings from their SM values are small. In 
Sec. IVIt we discuss the implications radiative corrections would have on our results. 

We finish in Sec. IVIIl by comparing the predictions of the individual models to each other 
and to the expected experimental sensitivity of the LHC and ILC. We present a decision 
tree for identifying the underlying model based on the couplings of h and point out which 
models cannot be distinguished even if the couplings of h were known exactly. Overall we 
find 15 models (or sets of models) with extended Higgs sectors that are distinguishable in 
principle in at least part of their parameter spaces. We also discuss the prospects for model 
differentiation based on the expected accuracy of Higgs coupling measurements at the LHC 
and ILC. For representative models, we plot the regions in which h would appear SM-like 
given the expected experimental uncertainties at the LHC (Fig. [5]) and ILC (Fig. [6]). We 
also provide a summary table showing the decoupling behavior of the h partial widths. We 
end with a brief summary of our conclusions. 

We find it most straightforward to organize the models on the basis of the structure of 
the Yukawa Lagrangian - in particular, the number of different Higgs doublet (s) involved 
in fermion mass generation. In Sec. IIIII we consider models in which the fermion masses are 
generated by the vev of only one Higgs doublet. Models with this characteristic are: 

• The SM, in which the Higgs sector consists of only one SU(2)i doublet that gives 
masses to the W and Z bosons and all the quarks and charged leptons. This is the 
simplest realization of the Higgs mechanism. 

• The SM extended with one or more singlet scalars [s], [so], 31, 32, 33|. This model 
yields an overall reduction in Higgs couplings due to doublet-singlet mixing and its 
phenomenology has been studied extensively. Models with similar light Higgs phe- 
nomenolo gy i nclude unparticle models in which Higgs-unparticle mixing can suppress 
couplings [M] and Randall- Sundrum models in which Higgs-radion mixing also leads 
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to reduced couplings [351. l36l. |37| 



The Type-I two Higgs doublet model (2HDM), in which only one doublet couples to 
fermions, while both doublets are involved in the generation of the W and Z boson 



masses 
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39| . The sharing of the vev between two doublets can have a dramatic 



effect on the couphng pattern of the Higgs boson. 

• A Type-I 2HDM extended with one or more singlet scalars. 

• A Type-I 2HDM extended with one or more additional doublets. 

In Sec. |IV]we consider models in which the fermion masses are generated by the vevs of 
two different Higgs doublets. Models with this characteristic are: 

• The Type-II 2HDM, in which one doublet generates the masses of the up-type quarks 
while a second doublet generates the masses of the down-type quarks and charged 



leptons 
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41 



42 



43 
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45| . This fermion coupling structure appears at tree level 



in the minimal supersymmetric standard model (MSSM), contributing to the great 
popularity of the Type-II 2HDM. In the MSSM, radiative corrections involving super- 
symmetric particles can induce potentially significant couplings of the bottom quarks 



to the "wrong" Higgs doublet 
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471]; while this feature formally puts the MSSM 



Higgs sector outside our requirement of natural flavor conservation, we nevertheless 
consider the features of this extension as well. 



A Type-II 2HDM extended with one or more singlet scalars 4^, |48|, |49|, |50|, l5l|, l52|, |53|, 



A Type-II 2HDM extended with one or more additional doublets 56|. 



The "flipped" and "lepton-specific" 2HDMs [621. 163l. I64l|. in which the coupling assign- 
ments of the two Higgs doublets to up-type quarks, down-type quarks, and charged 
leptons are varied relative to the usual Type-II 2HDM. In the flipped 2HDM, one dou- 
blet generates the masses of the up-type quarks and charged leptons while the second 
doublet generates the masses of the down-type quarks. In the lepton-specific 2HDM, 



Singlet extensions of the Higgs sector are also popular in supersymmetric models 



53 



55 



56 



57 



58, 



58 



6C 



61|. 
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one doublet generates the masses of both up-type and down-type quarks while the 
second doublet generates the masses of the leptons. We also consider extensions of 
these two models with additional doublets and/or singlet scalars. 

In Sec. |V]we consider models in which the fermion masses are generated by the vevs of 
three different Higgs doublets. Models with this characteristic are: 

• A "democratic" three Higgs doublet model (3HDM-D), in which one doublet generates 
the masses of up-type quarks, a second doublet generates the masses of down-type 
quarks, and a third doublet generates the masses of the charged leptons. 

• The 3HDM-D extended with one or more singlet scalars. 

• The 3HDM-D extended with one or more additional doublets. 

The physical Higgs boson content of these models can be summarized as follows. Consider 
a model that contains Nd complex doublets and A^^ singlets, Nc < Ns of which are complex. 
After removing the unphysical charged and neutral Goldstone bosons, this model contains 
Nd — l charged Higgs bosons H^, N^ + Ns CP-even neutral Higgs bosons H^, and Nd + Nc — l 
CP-odd neutral Higgs bosons A^. The state that we consider throughout is one of the CP- 
even neutral Higgs bosons, denoted h. Our results can be applied to any of the CP-even 
neutral Higgs bosons Hf. We make no assumptions about whether additional states can be 
discovered at the LHC. 

II. HIGGS COUPLINGS IN THE STANDARD MODEL AND BEYOND 

The Higgs doublet of the SM is given by 



where the vev of the Higgs field is vsm = 246 GeV. The couplings of the Higgs to SM 
fermions are given by the Yukawa Lagrangian, 




(1) 



^Yuk = -VeeR^^LL - yddR^^Qh - VuUr^^Ql + h.c. 



(2) 
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where Ql — (ul, di) , Ll — (i'l, cl) , and $ is the conjugate Higgs multiplet, 

$ = ^CT2$* =1 ^_ I . (3) 

These couphngs generate the fermion masses, ruf — UfVsM/V^, and the Feynman rules 
for the couphngs of the physical SM Higgs boson h — 0°''' to the fermions, — iy//-\/2 = 
-iruf/vsM = -iQf^- 

The couplings of the Higgs to the W and Z bosons arise from the covariant derivative 
terms in the Lagrangian, C = \T>^^\ , where the covariant derivative is given by 



V, = d^- t-^ + W-T-) - t./g^ + g'^Z^ (T^ - sin^ 0wQ) - leA^Q. (4) 

This term generates the W and Z boson masses, 



gvsM _ V9^ + 9 
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VSM 



mw = '^z = 2 ' 

and the Feynman rules for the couplings olhioW oi Z boson pairs, given by igv^g^v where 

SM _ g^^sM SM _ (g' + g'^>sM 

gw — 2 ' ~ 2 ' 

We now introduce our general framework for the extended models considered in this 
paper. In a model with multiple Higgs doublets, we can define the neutral, CP-even Higgs 
mass eigenstate under consideration as 

h = ^ ai(j)i, (7) 



where 0j = (f)^'^ is the properly normalized real neutral component of doublet ^j. The 
coefficients = {h\(f)i) are constrained by the usual quantum mechanical requirement that 
h be properly normalized: 

i 

In such a model, the vev that gives rise to the W and Z masses is in general shared among 
the doublets. We define the ratio of each doublet's vev to vsm as 

Vsm 
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where the second condition is required to obtain the correct W and Z masses. In the absence 
of CP violation, which we assume throughout, the quantities bi can all be chosen real and 
positive. Eq. [9]can also be thought of as a normalization requirement. In models containing 



only doublets, one can define a linear transformation to a "Higgs basis" [65|] in which only 
one doublet, carries a nonzero vev. The ratios hi are then given by bi = such 
that (py = bi(f)i and the condition Y2i bf = ^ follows from unitarity. 

In models containing both doublets and singlets, these relations are modified as follows. 
Because the Higgs state h can contain a singlet admixture, the sums in Eqs. [7] and [8] must 
include the singlet states as well as the doublets: 

h= ^ ai(j)i, ^ |ai|^ = 1, (10) 

doublets, singlets doublets, singlets 

where now (pi can also represent the real neutral component of a singlet. The W and Z 
masses, however, are generated only by the vevs of doublets, so that the sum in Eq. [9] is 
restricted to run over doublet vevs only: 

E b^ = ^- (11) 

doublets only 

While singlet scalars can have vevs of their own, these singlet vevs play no role in our analysis 
and will be ignored. 

In such an extended model, the couplings of the real neutral states (pi to W oi Z boson 
pairs arise from the covariant derivative terms for the doublets, 

C= Yl I^^AI'- (12) 

doublets 

After the mixing in Eq. (TDl the coupling of the physical state h to W ot Z boson pairs is 
controlled by the overlap of h with the doublet (p^ that carries the vev in the Higgs basis, 

9v = 9'v'{h\<P.), (13) 
where V = W oi Z. Inserting a complete set of states, we obtain, 

9v = = av"" E (14) 

i doublets only 

where the restriction of the sum to run over only the doublets arises because 0„ cannot 
contain a singlet admixture. Note that g^/dw^ = 9z/9z^^ ■ 
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We note here that the famihar 2HDM sum rule [69|] for Higgs couplings to gauge bosons 
can be generalized to models containing arbitrary numbers of doublets and singlets. Sum- 
ming over all CP-even neutral states and using Eq. [131 we have, 

= (^^n^E K^°l<^^')l' = (9^', (15) 

where the last equality is a consequence of completeness of the set of states H^. 



Under our assumption of natural flavor conservation [18|, ll9| , the masses of each type of 
fermion are generated by only one doublet. For fermion species /, the Yukawa Lagrangian 
can be written in the general form 

CYuk = -yffR^]FL+\i.C., (16) 

where Fi is the appropriate left-handed fermion doublet and $/ is the Higgs doublet that 
gives mass to fermion species /. If /i? is an up-type quark then $| should be replaced by <l>j 
above. These couplings generate the fermion masses, nif = yfVf/\^ = yjhfVsM/V^- Note 
that perturbativity of the Yukawa couplings, ?// ^ V^vr, together with the known third- 
generation fermion masses, imposes lower bounds on hf in each fermion sector: ^ 0.3, 
hd > 0.005, and hn > 0.003. After the mixing in Eq. [10], the coupling of the physical state h 
to an // pair is given by 

A sum rule can also be constructed for the fermion couplings, as follows: 

E(9;-V = fEiWWi^ = f = i^. (18) 

zrO rrO / SM 

i i 

where we have again used completeness of the set of states H^. Note that, unlike for the 
gauge coupling sum rule in Eq. [15], the right-hand side is not known a priori and instead 
depends on bf. We will show later that in some models, 6/ can be extracted from the 
couplings of a single state h. 

For compactness we define barred couplings normalized to their corresponding SM values, 

9^ - (19) 

so that (/j — > 1 in the SM limit and 

9w = ^ aA, 9.f = Y- (20) 

doublets ^ 
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III. FERMION MASSES FROM ONE DOUBLET 



In this section we consider extensions of the Standard Model in which the masses of all 
fermions arise from couplings to a single Higgs doublet. This class includes models containing 
additional doublets that do not couple to the fermions (so-called Type I models), as well as 
models containing one or more singlets. 



A. Standard Model plus one or more singlets (SM+S) 

The simplest way to extend the Standard Model Higgs sector is to add one real singlet 
scalar, S, which mixes with the usual SM Higgs boson to form two CP-even neutral Higgs 
mass eigenstates. The constraints of Eqs. [TOl [TT] become 

a} + al = 1, b} = 1, (21) 

where the subscript / refers to the Higgs doublet, which is responsible for all fermion masses, 
and s refers to the singlet. As noted below Eq. [H 6/ can be chosen real and positive, bf = 1. 
Simultaneously, aj can be chosen real and positive through an appropriate rephasing of the 
mass eigenstate h; can then be chosen real and positive through a rephasing of the field 
S. In particular, we can write 

af = v/r^ ^ v^. (22) 

The couplings of h to SM particles, normalized to their SM values as in Eq. [T9|, are then 
given by 

9w = ajbf = v^^, g^ = ^ = y/^. (23) 

In particular, the couplings of h to W or Z boson pairs and to fermion pairs are all scaled 
down by a common factor < 1 relative to their values in the SM. The production cross 
sections, partial decay widths, and total width of h are then all suppressed by a factor of ^ 
relative to those of the SM Higgs boson ^: 



SM 



e = l-<- (24) 



^ Note that a nonzero branching fraction of h to invisible particles could mimic this effect by suppressing 
Higgs event rates in all visible channels by a common factor. This possibility can be tested through a 



dedicated search for h invisible 



61 



66 



67 



68| 
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A measurement of any of these quantities allows a unique determination of the model pa- 
rameters a/ = and = y/1 — Because all the Higgs partial widths scale the same 
way with ^, the branching fractions of h are the same as in the SM. 

We also note that a decoupling limit can be defined in which — > and the couplings 
of h approach those of the SM Higgs; defining a small decoupling parameter 5 = we can 
write 

gf = gw = Vl^^l-l6\ ^ = 1-52. (25) 

These results can easily be extended to models containing two or more singlets by making 
the replacement 

E <■ (26) 

singlets 

The relations given above for af continue to hold, and we conclude that models with more 
than one singlet cannot be distinguished from the one-singlet model on the basis of the h 
couplings alone. 



B. Type-I Two Higgs Doublet Model (2HDM-I) 



The Type-I 2HDM |38l . l39l | has been extensively studied in the literature. This model 
contains two scalar SI]{2)l doublets, which we denote by $/ and $0; couples to fermions 
and $0 does not. 

The constraints of Eqs. [TOl [11] become 

a} + al = 1, b) + bl = 1. (27) 

The couplings of h to SM particles, normalized to their SM values as in Eq. [19], are then 
given by 

gw = cifbf + aobo, 9f = T~- (28) 

bf 

The vev ratios 6/ and 60 can both be chosen real and positive. Simultaneously, gw can be 
chosen real and positive through an appropriate rephasing of the mass eigenstate h. There 
is no freedom left to choose aj positive, though; thus gj can have either sign. 

Note in particular that while the couplings of h to fermions of all three sectors are scaled 
by the same factor relative to the SM, cju = (jd = ge, the coupling of h to W or Z boson 
pairs is scaled by a different factor unless bo = 0. This distinguishes the 2HDM-I from the 
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FIG. 1: Surface inhabited by the 2HDM-I in the plane of T^/F^^ versus Vl^/T^. Here tan/3 = 
Vf/vo = bf/bo. Note the double covering of the plane for small values of r^r-?^^ and V^/T^^ 
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SM plus a singlet discussed above. In the limit bo ^ we obtain gw = 9/ = cif and the 
couplings of h in the 2HDM-I reduce to those in the SM+S model. 

The constraint equations and coupling relations can be solved explicitly for the bi and 
factors in terms of the h couplings: 

1/2 



1 +6'/ - '^gwQf 
af = bfQf, 



ao 



9w - bjgf 



(29) 



Note that a full, unique solution is obtained if the relative signs oi cjf and cjw are known. If 
the relative signs are not known, then there are two possible solutions, as illustrated in Fig. [H 
Access to the relative signs of the couplings requires interfering the relevant amplitudes, e.g., 
by using h 'yy or h —>■ Zj. 

It is useful to make contact with the usual notation for the 2HDM-I 4J]. We assume 



that the state h under consideration is the lighter of the two neutral CP-even Higgs mass 
eigenstates, 

= \p2 (cos a Re $j — sin a Re <I>q) = cos a 0/ — sin a 0o, (30) 
with — 7r/2 < a < 7r/2, so that a/ = cosa and oq = — sin a. The ratio of the doublet vevs 
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can be defined according to 

tan/3 = ^ = ^, (31) 
^^0 bo 

so that bf = sin/3 and bo = cos (3. With these definitions, the couphngs of h to SM particles 
become 

cos (y 

Qw = sin(/5 — a), df = 7: = sin(/3 — a) + cot (3 cos(/5 — a). (32) 

sin p 

We note that for tan/5 < 1 the fermion couphngs can be enhanced {cjf > 1). Perturbativity 
of the top quark Yukawa couphng requires tan/5 > 0.3; there is no upper bound on tan/?. 
For fixed tan/5 the maximum possible value of F^/Fy*^ is 1 + cot"^ (3. 

The decoupling limit of this model occurs when the mass eigenstate h coincides with the 
state (pv, the vev-carrying doublet in the Higgs basis. In that case = {h\4>i) = {4>v\(pi) = h, 
so that (jf = af/bf = 1 and gw = cifbf + aobo = b'j + b^ = 1. Near the decoupling limit, we 
can parameterize the deviations of the couplings from their SM values in terms of a small 
parameter, 

S = cos(/5 — a) = ttfbo 
We have, for the couplings of h to gauge bosons, 

gw = Vl^^l-l6', 

The couplings of h to fermions depend also on tan/5: 

(jf = Vl - 52 + cot /5 5 ~ 1 + cot /5 5 - 

where the terms of order must be kept if cot /5 is very close to zero. Note that 8 can take 
either sign. 



aobf. (33) 



J- 



J- w 



^ (34) 



~ l + 2cot/5(5-(52, (35) 



C. 2HDM-I plus one or more singlets (2HDM-I+S) 

We next consider the consequences of adding a real singlet scalar field, S", to the 2HDM-F 
The constraints of Eqs. [TOl [11] become 

aj + + = 1, b)^bl = 1, (36) 
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where = {h\S) and the other Cj, 6j are defined as in the previous section. The couphngs 
of h to SM particles, normahzed to their SM values, are given as for the 2HDM-I by 

a f 

gw^CLfbf + aobo, 9f = t- (37) 

As before, bf, bo, and gw can be chosen real and positive, while gf can have either sign. The 
coefficient can then be chosen real and positive by a rephasing of S. 

With five parameters and only four equations, the parameters of this model cannot be 
fully solved for in terms of the h couplings. In order to display the ambiguity we define 

C = l-al = a} + al, (38) 

with < ^ < 1 parameterizing the doublet content of h. We then obtain, 

1/2 



^ + 9f- '^9wgf 



bo-Jl-bj, 



9w - bjgf 



«/ = bfgf, ao = — ^ as = a/I -C, (39) 



1-b' 



where ^ remains an undetermined parameter. 

This model can be cast into the usual notation of the 2HDM-I as follows. We first 
parameterize the doublet-singlet mixing in terms of ^, 

h^^/^h' + ,/l^S, (40) 

where h' corresponds to our Higgs in the 2HDM-I in the limit of zero singlet admixture: 

/i' = cosq;0/ — sina^o- (41) 

We then have af — cos a, ao — sin a, and as — ^/l — ^. The couplings are given by 

/7 V /-COS a 

9w = sm(/3 -a), 9f ^ — ^- (42) 

In particular, the couplings of h to SM particles are all scaled down by a common factor V^- 
Regardless of whether the parameters of the model can be solved for uniquely in terms 
of the h couphngs, this model would be distinguishable from the 2HDM-I if it occupied 
a different footprint in the space of observables; i.e., if one could obtain sets of couplings 
{gw,9f) in this model that could not be obtained in the 2HDM-I. This is not the case. Any 
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set of couplings {gw,gf) that can be obtained in the 2HDM-I+S can also be obtained in 
the 2HDM-I, albeit from different underlying values of the parameters a/, ao, and bo. In 
particular, the models are identical when = 1; away from this limit the ellipses in Fig. [1] 
are simply scaled down by ^ on both axes. 

The presence of the singlet thus cannot be established through measurements of the 
couplings of h only. However, if the couplings of a second CP-even neutral Higgs state (H) 
could be measured, nonzero singlet mixing would violate the usual 2HDM coupling sum 
rule [69|], (g^)^ + {gwY = idw^y (see Eq. [15]). This violation would indicate the presence 
of a third CP-even state such that Y^Liiawf = iaw^Y- 

These results can easily be extended to models containing two or more singlets by making 
the replacement 

E < = l-^- (43) 

singlets 

Again, such a model cannot be distinguished from the 2HDM-I on the basis of the h couplings 
alone. 

We conclude that adding one or more singlets to the 2HDM-I results in a model that 
cannot be distinguished from the 2HDM-I on the basis of the h couplings alone. 



D. 2HDM-I plus additional doublet(s) (2HDM-I-HD) 

Let us now consider the consequences of adding one or more additional Higgs doublets to 
the 2HDM-I. These additional doublets can carry vevs but, under our assumption of natural 
flavor conservation, they must not couple to fermions. We can denote the field content of 
the model as $/, ^oi, with i = 1 . . .n {n > 2) counting the doublets that do not couple to 
fermions. 

We first define a linear combination 0q of the neutral CP-even states such that 

h = af(j)f + ao0o, aj + = 1. (44) 

The vev of (pQ is parameterized by 6q = {(p^lcpy), chosen to be positive; the phase of h is 
chosen to make gw = (^fbf + o'o^'o positive. Eq. [11] is modified to read 

6j + b'^ =uj^, 0<uj< 1, (45) 
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where u < 1 indicates that a nonzero vev is carried by the hnear combination(s) of (/)0i 
orthogonal to h. Again, this model has five parameters but only four constraint equations; 
the solution for the model parameters becomes 

1 1/2 



5/ 



1 + u'^gf - 2gwgf 
bf9f, a'o = 



9w - b}gf 



(46) 



where uj remains an undetermined parameter. 

Translating into the usual 2HDM-I notation, we define a as for the 2HDM-I and tan j3 as 



n bf 
tan/3 = -f , 

where now sin/3 = hf/uo and cos/3 = Vq/uj. In this notation, the couplings of h become 

1 cos a 



gw = ^ sin(/3 — a). 



9} 



(47) 



(48) 



u sin (3 

The couplings oihto gauge bosons are scaled down by a factor u; < 1 while the couplings of h 
to fermions are scaled up by a factor l/oj > 1. Again, this model occupies the same footprint 
in Higgs coupling space as the 2HDM-I (consider Fig. [1]), and we conclude that adding one 
or more additional doublets to the 2HDM-I results in a model that cannot be distinguished 
from the 2HDM-I on the basis of the h couplings alone. This conclusion remains unchanged 
if singlets are added to the model as well. 



IV. FERMION MASSES FROM TWO DOUBLETS 



We now consider models in which the fermion masses arise from couplings to two different 



Higgs doublets. Imposing natural fiavor conservation a 



couplings of two Higgs doublets to the fermions 



6d, 



63|, 



lows for three possible patterns of 



6|: 



(i) the Type-II 2HDM, or Model II, in which one doublet generates the masses of the 
up-type quarks while the other generates the masses of the down-type quarks and 
charged leptons (this is the coupling structure present at tree level in the MSSM); 

(ii) the "fiipped" 2HDM, in which one doublet generates the masses of the up-type quarks 
and charged leptons while the other generates the masses of the down-type quarks; 
and 
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(iii) the "lepton-specific" 2HDM, in which one doublet generates the masses of up- and 
down-type quarks while the other generates the masses of the charged leptons. 

We consider here also extensions of these three models obtained by adding one or more 
electroweak singlets or doublets that do not couple to fermions. 

While the Higgs sector of the MSSM is a Type-II 2HDM at tree level, one-loop ra- 
diative corrections involving supersymmetric particles can induce a significant coupling of 
the bottom quark to the "wrong" Higgs doublet, encoded in an extra coupling parameter 



7l\ . This violates our assumption of natural flavor conservation; however. 



for completeness, we consider the main features of this model here separately. 



A. Type-II Two Higgs Doublet Model (2HDM-II) 

The Type-II 2HDM Q, B, Q is perhaps the most widely studied extension of the SM 
Higgs sector. The Higgs content and coupling structure are the same as in the MSSM at tree 
level. This model contains two scalar SU(2)i doublets, which we denote by and ^d- 
generates the masses of the up-type quarks while $rf generates the masses of the down-type 
quarks and the charged leptons. 

The constraints of Eqs. [10] and [TT] become 

al + al = l, bl + bl = l. (49) 

The normalized couplings of h to SM particles are then given by 

gw = ajju + adbd, Qu = 7^, 9d = 9e = 7^, (50) 

bu bd 

where (ju, Qd, and gi denote the normalized couplings of h to all three generations of up-type 
quarks, down-type quarks, and charged leptons, respectively. 

The vev ratios b^ and b^ can both be chosen real and positive. Simultaneously, we can 
choose gw to be real and positive through an appropriate rephasing of the mass eigenstate 
h. There is no freedom left to choose the signs of the fermion couplings (ju and gd', depending 
on the underlying values of the parameters they can take the signs ++, H — , or — h; gu and 
(jd cannot both be negative. 
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The 2HDM 
relation I 



-II has four parameters related by five constraints, resulting in a pattern 
751 1 among the three couplings oi h:^ 



ud 



9w{9u + 9d) ~ 9u9d = 1- 



(51) 



This pattern relation provides a test of the 2HDM-II coupling structure: It defines a two- 
dimensional surface accessible by the model in the three-dimensional space of couplings gw, 
gu, and (jd- 

The constraint equations and coupling relations can be solved explicitly for the 6, factors 
in terms of the h couplings: 

1/2 



bu 
b, 



gw 


- 9d 


1/2 


" 1 


-9d ' 




-9d_ 




.91 


-9'd\ 


'gw 


- 9u 


1/2 


' 1 






- gu _ 




[fd 


-9l. 



(52) 



where in the second relation the dependence on g^ has been removed using the pattern 
relation. We also obtain the a,- factors, 



bugui 



ad = bdgd- 



(53) 



If the relative signs of (ju, (jd, and gw are known, then the solution for the model parameters is 
unique. Note also that, unlike in the 2HDM-I, a unique solution for bu and bd can be obtained 
even if the signs of the couplings are not known, by using the second set of equalities in 
Eq. [321 In the absence of information on the signs of the couplings, the magnitudes of au 
and ad can also be determined uniquely but their relative signs cannot. 

Let us now make contact with the usual notation for the 2HDM-II 4J]. We assume 



that the state h under consideration is the lighter of the two neutral CP-even Higgs mass 
eigenstates, 

/i° = v^(cosaRe$° -sinaRe$S), (54) 

so that au = cos a and = — sin a. The ratio of the doublet vevs can be defined according 
to tan/? = Vu/vd = bu/bd, so that bu = sin/3 and bd = cos/?.'' With these definitions, the 



^ One can define an additional pattern relation involving gw, gu and gi according to Put = gw{gu + gi) — 
9uge = 1 = Pud- 

^ The constraint gw > corresponds to (3 — w < a < f3. Perturbativity of the top quark Yukawa coupling 
requires tan /3 > 0.3; perturbativity of the bottom quark Yukawa coupling requires tan (3 < 200. 
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FIG. 2: Surface inhabited by the 2HDM-II in the plane of (left) T^/rf^^ and (right) T^/T^'^ 
versus r(y/r^, for various values of tan/?. 

couplings of h to SM particles become 

gw = sin(/?-a), 
cos a 



9u 
Qd 



sin (3 



sin(/3 — a) + cot (5 cos(/3 — a) 



sm a 

cos [3 



sin(/5 — a) — tan [3 cos(/5 — a). 



(55) 



We note that tan (3 can be obtained from coupling measurements using 

tan/3 = fil^ = . (56) 

The relations among the Higgs partial widths into up-type quark, down-type quark (or 
charged lepton), and W boson pair final states are shown in Figs. [2] and [3] for various values 
of tan/5. The key difference between the 2HDM-II and the 2HDM-I is the different behavior 
of cju compared to cjd (and gp), as illustrated in Fig. [31 In particular, the 2HDM-I would fall 
on a line of slope +1 through the SM point (1, 1) on this plot. 

The decoupling limit of this model occurs when the mass eigenstate h coincides with the 
state 0^,, the vev-carrying doublet in the Higgs basis. Near the decoupling limit we param- 
eterize the deviations of the couplings from their SM values in terms of a small parameter 



6 = cos(/5 -a) = ayhd - ajhu- 



(57) 
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FIG. 3: Surface inhabited by the 2HDM-II in the plane of T^/rf*^ versus T^/T^ 



itSM 



The couplings and corresponding partial widths, normalized to their SM values, become 

l-5\ 



gw = ~ 1 - 



J- w 



2" ' 



= Vl - 52 + cot /? (5 ~ 1 + cot /5 5, 



9d = de = Vl — S^ — tan (5 ~ 1 — tan /? (5, 
Note that 5 can take either sign. 



rSM 
J- w 



u 

-*- d 
p5M 



1 + 2cot/?5, 
l-2tan/?(5. (58) 



B. 2HDM-II plus one or more singlets (2HDM-II+S) 

We now consider the consequences of adding a real singlet scalar field, S", to the 2HDM-II. 
The constraints of Eqs. [TOl [11] become 



2 I 2 I 2 -1 
+ «d + «s = 1, 



(59) 



where = (/il^) and the other Oj, fcj are defined as in Sec. IIV A[ The couplings of h to SM 
particles, normalized to their SM values, are given in terms of and fo^.d as for the 2HDM- 
II by Eq. [50l As before, hd-, and gw can be chosen real and positive. The coefficient 
can be chosen real and positive by a rephasing of S. 

Because of the presence of the additional parameter a^, the pattern relation of the 2HDM- 



20 



II no longer holds. Instead we obtain 



Pud = 9w{9u + 9d) - guQd = ^ < 1, 



(60) 



where ^ = 1 



+ parameterizes the doublet content of h. In particular, ^ can be 



determined by applying the pattern relation to measurements of the couplings gw, Qu and 
Qd- The solutions for the rest of the model parameters then become 



K 
bd 



9w 


- 9d 


1/2 




-fd] 




-9d_ 




.91 


-9l\ 


gw 


-Qu 


1/2 






Jd 


- 9u _ 






- 9^ 



1/2 



1/2 



bu9u'i 



(I'd — bd9d, 



(61) 



where in the expressions for b^, bd we have shown how the dependence on can be traded 
for dependence on ^ using Eq. [601 

Clearly, if the relative signs of gw, (ju and cjd are known, then this model can be distin- 
guished from the 2HDM-II using the pattern relation {Pud is equal to one in the 2HDM-II 
and less than one in the 2HDM-II+S) and the solution for the model parameters is unique. 
If, however, the relative signs of gw, 9u and gd are not known, then ^ cannot be obtained 
uniquely and there will be discrete ambiguities in the solutions for all the parameters. In 
this situation the pattern relation can still be used to test for the presence of the singlet in 
the model; if no combination of signs of the Higgs couplings gives Pud = 1, then the model 
cannot be the minimal 2HDM-II. 

This model can be cast into the usual notation for the 2HDM-II as follows. We first 
parameterize the doublet- singlet mixing in terms of ^, 



(62) 



where h' corresponds to the Higgs state considered in the 2HDM-II in the limit of zero 
singlet admixture: 

h' = cos a(j)u — sin a (j)d. (63) 

We then have au = cos a and ad = — a/^ sin a. The couplings are given by 

■ cos a A- sin a 



9w= sin(/3 - 



a) 



9 



sin [3 



9d = 9e 



-VI 



cos (3 



(64) 
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In particular, the couplings of h to SM particles are all scaled down by a common factor 
< 1. This means that the 2HDM-II+S lives on a volume in the three-dimensional 
parameter space of gw^ (ju, and cjd, consisting of the surface inhabited by the 2HDM-II 
(corresponding to ,^ = 1) together with all lines that connect points on that surface to the 
origin (corresponding to < ^ < 1). Clearly, the 2HDM-II-I-S occupies a different footprint 
in couphng space than the 2HDM-II, and it can thus be distinguished from the 2HDM-II. 
This is different from the case of the 2HDM-I-I-S; the reason is that the Type-II fermion 
coupling structure yields a third observable coupling related nontrivially to the other two. 

The decoupling limit comprises 5 = cos(/5 — a) = a,J)d — ad^u ^ and e = y^l — ^ — > 0. 
(Note that S can have either sign while e is chosen positive.) The couplings and corresponding 
partial widths, normalized to their SM values, become 

i w 

, , rh 

= [y/T^ + cotp6] ^r^~ 1 + cot/35- ie^, ~ 1 + 2cot/35-e^ 

u 

g^^ g^^ [v^n^-tan/35] v^r^~ 1 -tan/35- ie2, ^ ~ i _ 2 tan/3 5 - e^. 

Note that in the hmit 5 — > with e finite, the deviations of the h couplings from their SM 
values become identical to those in the SM+S. 

These results can easily be extended to models containing two or more singlets by making 
the replacement 

singlets 

Such a model cannot be distinguished from the 2HDM-II-I-S (with only one singlet) on the 
basis of the h couplings alone. 

C. 2HDM-II plus additional doublet (s) (2HDM-II+D) 

We now consider the consequences of adding an additional Higgs doublet $o to the 
2HDM-II. The additional doublet can carry a vev, but under our assumption of natural 
flavor conservation it must not couple to fermions. The constraint equations become, 

al + al + al = l, bl + bl + bl = l, (67) 
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where Oq = (/i|0o) and Bq = Vq/vsm- The normahzed couphngs of h to SM particles are 
given by 

gw = ciubu + CLdhd + flo^o, (ju = 7^, 9d = ge = (^8) 

All three 6j parameters can be chosen real and positive; gyy can also be chosen positive 
through an appropriate rephasing of h. Any combination of signs is then possible for cju and 
Qd] in particular, both can be negative (for ao&o > l^ubu + ddbdl) in contrast to the 2HDM-II 
or 2HDM-II+S. 

Because of the presence of the two additional parameters oq and Bq, the model is undercon- 
strained and the parameters a,, bi cannot be extracted in terms of the h couplings. However, 
the model is distinguishable from the 2HDM-II because the pattern relation of Eq. [51] no 
longer holds. In some parts of parameter space, this model can also be distinguished from 
the 2HDM-II+S. 

In order to illustrate these features, we cast the model into the usual notation for the 
2HDM-II. We first parameterize the mixing of the third doublet in terms of an angle 6, 

h = cos9 h' + sinO (j)o, (69) 

where h' = cos a 0„ — sin a (pd corresponds to the Higgs in the 2HDM-II in the limit of 
zero mixing with the extra doublet. We then have a„ = cos cos a, = — cos 6' sin a, and 



ttQ = sinO. We also define tan/5 = Vu/vd = bu/bd and cosil = \/b\ + b\, sini7 = 60, where 
the angle < i7 < 7r/2 parameterizes the amount of vev carried by $0- The couplings of h 
are then given by 

= cosi7 COS0 sin(/3 — a) + sinfi sin^, 

cos6' cos a _ _ cos6' sin a , ^ 

9u = FT^IT' 9d = ge = ^- (70) 

cos i 2 sm p cos i 2 cos p 

We note the features of two limiting cases: 

(i) When bo = (i.e., cosfl = 1), the h couplings reduce to those of the 2IIDM-II+S, 
with a/^ replaced by cos 6*. This happens because in this limit, 0o does not couple to 
fermion pairs or gauge boson pairs and the physics is simply that of the 2HDM-II with 
mixing of a "sterile" state into h. The pattern relation in this special case becomes 

Pud = gw{9u + Qd) - guQd = cos^ 9 <l. (71) 
(Note that the 2HDM-II result Pud = 1 is recovered in the limit cos^ 6^1.) 
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(ii) When ao — (i.e., cos^ — 1, so h — h') there is no mixing of the new doublet into 

h, but the vev of 4>o is nonzero so that the total vev carried by the two doublets that 

couple to fermions is reduced. The fermion Yukawa couplings must thus be enhanced 

in order to yield the required fermion masses, while the coupling oi h to W or Z pairs 

is suppressed. In this case the couplings of h become 

_ 1 cosa _ 1 sin a 

gw = cos il sm{/3 - a), ^„ = n^^^ 9d^ p: ^, (72) 

cosiZ sm/3 cosU cosp 

and the pattern relation in this special case becomes 

sm cy, cos ex 

Pud = gw{9u + 9d) - QuQd = 1 + tan^ ^^-p, ^- (73) 

sm /? cos p 

In particular, Pud > 1 whenever sin a cos a > 0, i.e., whenever ga and g^ have opposite 
signs. Furthermore, when sin a cos a < (i.e., when g^ and cju have the same sign), 
small values of cosO, sin (3, and/or cos/? can yield Pud < 0. Either of these situations 
allows the 2HDM-II+D to be distinguished from both the 2HDM-II and the 2HDM- 
II+S. (Note that the 2HDM-II result Pud = 1 is recovered in the hmit cosQ — > 1.) 

In the general case of both cos 6 < 1 (nonzero mixing of 0o into h) and cos Q < 1 (nonzero 
vev of 0o) , the footprint of the 2HDM-II+D covers a three-dimensional volume in the space 
of couphngs {gw: gu: Qd)- Part of this volume overhes the footprint of the 2HDM-II-I-S (when 
mixing dominates, or when gu and ga have the same sign), but part is unique to the 2HDM- 
II+D (when vev sharing dominates, or when g^ and gd are both negative as discussed before). 
Thus the model is distinguishable in general from the 2HDM-II, and is distinguishable from 
the 2HDM-II+S in some regions of parameter space. 

We now describe the approach to decoupling in this model. The decoupling limit corre- 
sponds to {h\(l)y) — > 1. Deviations from this limit can be parameterized by writing 

h = c\\4>y + c±4>±, {4>±\4>v) = 0, (74) 

where c^^+ c]_ — 1 and (f)y is given in our parameterization by 

0^ = cos O, (sin P(f)u + cos P (f)d) + sin Q 0o- (75) 

The component of h orthogonal to (f)^ can be constructed as follows. We first define two 
states orthogonal to 0„ and to each other: 

(j)±i = cos /3 0u - sin P (76) 
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which hes in the (l)u^(pd plane, and 

cf)±2 = — sin f2 (sin P 4>u + cos /? (pd) + cos Q 0o- (77) 
Then can be parameterized in terms of a new mixing angle 7, 

0_L = sin 7 0_L1 + cos 7 0_L2 

= (sin 7 cos P — cos 7 sin Q sin /3) 0„ + (— sin 7 sin /3 — cos 7 sin f2 cos P) (pd 

+ (cos7 cosi7) 00- (78) 
Defining the decoupling parameter 5 = c_l, we obtain the couplings of h: 



gw = {h\(t)^) = Vl - 52 
K 



9d = 9i = —r = VI - 0^ + 6 



■ . cot/? ■ 

sm7 — ■ — cos 7 tan iZ 

cos \l 

. tan/3 ■ 

— sm 7 — ■ — cos 7 tan i 2 

cos 



(79) 



Letting 5 take either sign, we can fix < 7 < vr. Note that for sin 7 = cosfi = 1, these 
formulas reduce to those for the 2HDM-II given in Eq. [581 The decoupling limit corresponds 
to 5 ^ 0. 

This analysis can be extended to the 2HDM-II plus two or more doublets in a straight- 
forward way. We denote the doublets that do not couple to fermions as $oj, with i = 1 . . .n 
{n > 2). As in Sec. IIIIDl we first define a linear combination </)q of the neutral CP-even 
states (po- such that 

h = au4>u + ad(j)d + ao0o, + + = 1. (80) 

The vev of 0q is parameterized by 6q = {4>'q\4>v), and Eq. [11] becomes bl^ + bd + b'^ < 1, where 
the inequality accounts for the vev carried by the linear combination(s) of (pQ. orthogonal 
to h. While the underlying physics of this model differs from that of the 2HDM-II plus one 
extra doublet, the footprint of the model in coupling space is the same. This can be seen 
straightforwardly by noting that the pattern relation P^d can take any value in the 2HDM- 
II+D, leaving no room for a larger footprint when additional extra doublets are added. 
Thus, while it is possible to tell from the couplings of h alone that (at least) one additional 
doublet has been added to the 2HDM-II, it is not possible to tell how many. 

The addition of singlet (s) to the 2HDM-II+D can be parameterized in a similar way. We 
first note that as far as the couplings of h are concerned, adding an additional Higgs doublet 
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with zero vev is indistinguishable from adding a singlet. We again obtain Eq. [HO] in which 
0Q now denotes the appropriate linear combination of (po and the singlets. The constraint 
equation for the doublet vevs remains as given in Eq. [671 We see that the model occupies 
the same footprint in h coupling space as the 2HDM-II+D and thus it is not possible on the 
basis of h couplings alone to tell whether the 2HDM-II+D also contains additional singlets. 



D. Flipped 2HDM, lepton-specific 2HDM, and their extensions 



The flipped and lepton-specific two Higgs doublet models comprise the two possible al- 



ternate assignments of fermion coup 



models were introduced in Refs. 
have been made in Refs. 



76 



62 



i ngs of the two-doublet models considered here. These 
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64| . Some early studies of their phenomenology 



77( 1 . Much can be extrapolated in a straightforward way from 



existing results for the usual 2HDM-I and 2HDM-II. 

In the flipped 2HDM, one doublet generates the masses of the up-type quarks and the 
charged leptons while the other doublet $rf generates the masses of the down-type quarks. 
The constraint equations remain identical to those of the 2HDM-II as given in Eq. HHl while 
the normalized couplings of h to SM particles are given by 



gw — o,uhu + adhd, 



a^ 



ad 

9d = -7-- 
bd 



(81) 



The distinguishing feature of this model is the behavior of gi. The quark coupling results 
carry over unchanged from the 2HDM-II model and its extensions by extra doublets and/or 
singlets. 

In the lepton-specific 2HDM*, one doublet $g generates the masses of all flavors of quarks 
while the other doublet $f generates the masses of the charged leptons. The constraint 
equations become 



2 2 

a„ + ag 



hi + 



1. 



The normalized couplings of h to SM particles are given by 

gw = aqbq + aibi, gu = gd = 



9i 



ae_ 



(82) 



(83) 



® LHC phenomenology for h in the lepton-specific 2HDM was also discussed in Ref. 7^. Ref. [79| also 
makes use of this fermion coupling structure. 
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Note that in the quark sector, this model is identical to the 2HDM-I. Its distinguishing 
feature, however, is again the behavior of gf, the pattern relation and all other results for 
the 2HDM-II and its extensions by extra doublets and/or singlets carry over to this model 
with the replacements 

9u 9q, 9d 9i- (84) 



E. MSSM (2HDM-II with A^) 



At tree level, the Higgs sector of the MSSM is a Type-II 2HDM. The natural flavor con- 
servation structure of the Yukawa couplings is enforced by the analyticity of the superpo- 
tential. Beyond tree level, however, radiative corrections involving loops of supersymmetric 
particles can induce additional couplings of right-handed fermions to the "wrong" Higgs 
doublet 4^, 47, 70, 71 1. Thus, beyond tree level the Higgs sector of the MSSM is techni- 
cally a Type HI 2HDM^. This violation of natural flavor conservation is a consequence of 
supersymmetry breaking; because of this, the loop-induced wrong- Higgs couplings do not 



82 



83|. 



decouple as all SUSY mass parameters are simultaneously taken large 

The most important loop-induced wrong Higgs couplings of this type arise in the bottom- 
quark sector from loops involving bottom squarks and gluinos (involving the large QCD 
gauge coupling) and from loops involving top squarks and charginos (involving the large top 
Yukawa coupling). In particular, the effect of the wrong-Higgs coupling on gi, is enhanced by 
tan/?, meaning that even though it is a one-loop effect, it can be important at large tan (3. 
The radiatively-corrected couplings can be parameterized by an effective Lagrangian 8^] , 



- C,s = e,,hbnH',Qi + AhbRQlH^* + h.c 



(85) 



with Hu and Hd defined in the usual way for the MSSM with opposite hypercharges. Note 
that we absorb into hh any SUSY radiative corrections to the "right-Higgs" coupling. 
The physical bottom quark mass is given by 



rUb 



V2 V2 



hbVsM cos P 
V2 



1 + 



Ahbtan(3\ _ hbVsM cos P 



hh 



V2 



Ab). 



(86) 



^ The phenomenology of the general Type III 2HDM has been reviewed in Ref . |80| . In this model the basis 
chosen for the two Higgs doublets is somewhat arbitrary; basis-independent methods have been developed 



Refs. 
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(Note the factor of tan/3 that is absorbed into the definition of A;,.) Similarly, the h^bb 
coupling becomes 



9^ 



sm a 



h 

where the mixing angle a is defined as in Eq. 
written in terms of rrib and A;, by noting that 

nib 1 



cosa- 



Ahb 
V2 



(87) 



for the 2HDM-II. This coupling can be 



h 

VSM COS 1 + Afe ' 



nib 



A, 



x/2 vsM sin ^ 1 + Afe 



Inserting these relations into and normalizing by the SM coupling we obtain 



gb = sin(/? — a) — tan P cos(/? 



a] 



1 - cot^ /3 Ab 



(88) 



(89) 



l + Ab 

The Ab corrections are typically the only large SUSY radiative corrections to the Higgs 



Yukawa couplings 



84| - in particular, the analogous corrections to the Higgs couplings 



to top quarks are not tan/3 enhanced, and those to the Higgs couplings to tau leptons 



involve only the small electroweak gauge couplings 



7l|. Thus the SUSY corrections to these 



couplings can generally be neglected, and the usual 2HDM-H relations are recovered: 



gw = sin(/?-a), 
cos a 

gu 
9e 



sin P 
sin a 

cos P 



sm{P — a) + cot P cos(/? — a), 
= sm{P — a) — tan P cos(/5 — a). 



(90) 



The model parameters can be obtained as in the 2HDM-II by using the couplings that 
are unaffected by A;,: 



tan P 



7,2 



- 1 



nl/2 



l-9l 



cos a = gu 



9l - 9} 



sm a = gi 



84| 



1 -2 -\ 1/2 
^-9V 



fi-9l 



(91) 



(92) 



The value of A^ can also be extracted from the h couplings using 

. 9h-gi 
IXb — ^• 

Qu - Qb 

We first note that the pattern relation of the 2HDM-II among the couplings gw, gu and 
ge survives: 

Pu£ = 9w{gu + gt) - 9u9t = 1- (93) 

This allows the MSSM Higgs sector to be distinguished from the more general three-Higgs- 
doublet models discussed in the next section and allows one to test for the presence of 
additional singlets or doublets that mix with h or carry nonzero vevs. 
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However, the A;, correction to the Higgs couphng to bottom quarks leads to ^ ge, such 
that the pattern relation among the W, u and d couplings is violated: 

, Afefl + COt^/9) , ^, 

Pud = gw[9u + gd) - gugd = l - cos [(3- a) — — . (94) 

Depending on the sign of A;,, the right-hand side can be greater or less than one. 

In the decoupling limit the deviations of the h couplings from their SM values can be 
parameterized in terms of 5 = cos(/3 — a). We have, 

J- w 



g^ = ^l- 5"^ + cot 135 ~ 1 + cot/55, ~ 1 + 2cot^(5. 



u 



gb = - 1 - tan/?' 5, - 1 - 2tan/5'(5, 



= v/r3^_tan/55 ~ 1 -tan/35, ^ 1 - 2tan/5(5. (95) 



Note that 5 can take either sign. These expressions are identical to those given in Eq. 
for the 2HDM-II except that in the bottom quark couplings we have replaced tan (3 with 

tan/3' = tan/3 1"^"^^^^ (96) 
1 + Ab ^ ^ 

V. FERMION MASSES FROM THREE DOUBLETS 

Finally we consider a model in which the fermion masses arise "democratically" from 
couplings to three different Higgs doublets - i.e., models in which the masses of the up- type 
quarks, down-type quarks, and charged leptons are generated by couplings to three different 



Higgs doublets ^d, and $£, respectively. Such a model was considered in Ref. j6J]. We 
also consider extensions of this model obtained by adding one or more singlets or doublets 
that do not couple to fermions.^° 

A similar Higgs-fermion coupling structure has recently been proposed in the "Private 



Higgs" model 



85 



861], which introduces one Higgs doublet for each of the six flavors of 



quarks in order to address the hierarchy of quark masses. Here, however, we limit the 



Models with three or more doublets introduce the possibility of new CP-violating parameters in the n x n, 
n > 3 mixing matrices of the charged scalars [3] ■ Again, we neglect the possibility of CP-violating effects 
in this work. 
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discussion to models in which the masses of fermions of a given electric charge are generated 
by their couplings to only one Higgs doublet; in particular, we do not allow the Higgs 
coupling structure to differ by fermion generation. We also make no assumptions about the 
structure of the Higgs potential. The discussion here can be extended to models in which 
the three generations are treated differently but care must be taken to avoid Higgs-mediated 
flavor- changing neutral currents. 



A. Democratic Three Higgs Doublet Model (3HDM-D) 



In this model the up-type quarks, down-type quarks, and charged leptons each get their 
mass from a different Higgs doublet, denoted $„, and respectively. The constraints 
of Eqs. [TOl [11] become 



o-l + + a« = 1, 



bl + bl + b^, = l. 

The normalized couplings of h to SM particles are then given by 

gw = ciubu + adbd + a-ibi, 



9u = 7-, 
bu 



ad 
bd 



(97) 



9i = ^. (98) 



The vev ratios bu, bd, and be can all be chosen real and positive. Simultaneously, we can 
choose gw positive through an appropriate rephasing of the mass eigenstate h. There is no 
freedom left to choose the signs of the fermion couplings (ju, Qd, or gi, depending on the 
underlying values of the parameters they can take any combination of signs so long as at 
least one of them is positive. 

The constraint equations and coupling relations can be solved explicitly for the b-i factors 
in terms of the h couplings: 



bu = 


'1 - gwigd + Qi) 


+ gdge' 


1/2 


. {gu-gd){gu- 


-gd . 




bd = 


1 - gw{gu + gi) 


+ guge 


1/2 


. {gd-gu){gd- 


-gd . 




bi = 


1 -gw{gu + gd] 


+ gugd 


11/2 


. {gt - gu) [ge - 


-gd) . 





(99) 



We also obtain solutions for the factors, 

au = buCju, ad = bdcjd. 



ae = bege. 



(100) 
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As in the other solvable models, if the relative signs of g^, Qdi Oe and cjw are known, then the 
solution for the model parameters is unique. However, if the relative signs are not known, 
discrete ambiguities arise in the solutions for the 6j and factors. 

The key feature that distinguishes the democratic 3HDM from the previous models con- 
sidered is that Qu ^ Qd ^ ge- While this is also true in the MSSM with Af, corrections, the 
MSSM couplings satisfy the pattern relation P^i = 1 involving the couplings gw, du, and gi 
(Eq. [93]); this relation does not hold in the 3HDM-D. 

The features of this model can be clarified by examining its parallels with the 2HDM- 
II+D. In particular, the behavior of cjw^ 9u and cjd is identical to that in the 2HDM-II+D, 
while now gg behaves differently with the model parameters and serves to distinguish the 
current model. As in the 2HDM-II+D, we parameterize the mixing according to h = 
cos Oh' + sin 9 (j)£, with h' = cos a(f)u — sin a 0^, yielding = cos 6 cos a, = — cos 6 sin a, 
and ai = sin 6*. We also define tan/5 = Vu/vd = bu/bd and cosi7 = + fc^, sinf2 = be. The 
couplings of h are then given by 

gw = cosficos6'sin(/9 — a) + sinnsin6', 

cos^ cos a _ cos^ sin a _ sin^ 

9u = n—^^ 9d = p: ^, 91 = (101) 

cos \ I sm p cos i I cos p sm \ I 



The decoupling relations can be parameterized exactly as in the 2HDM-II+D (Eq. [79 
except for gi, which is given by 



9e 



Vl- 52 + [cos 7 cot Q], (102) 



bi 

where (pe is the third doublet and 7 is defined as in Eq. [78] with 



B. 3HDM-D plus one or more singlets (SHDM-D+S) 

We now consider the consequences of adding a real singlet scalar field, S, to the 3HDM-D. 
The constraints of Eqs. [TOl [11] become 

al + al + a\ + = 1, + &rf + b\ = 1. (103) 

The formulae for the normalized couplings of h to SM particles are identical to those of the 
3HDM-D given in Eq. [98] The vev ratios bu, bd and be and the coupling gw can all be chosen 
real and positive; can then be chosen real and positive by an appropriate rephasing of S. 
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Because of the presence of the additional parameter a^, this model is distinguishable from 
the 3HDM-D in part of its parameter space, as we now show. First we define the following 
three combinations of h couplings, 



Xu = 
Xd = 
X, = 



gw{9d + ge) +gdge 



{(ju - gd){gu - 


-gi) 


1 - gw{gu + gd 


+ gug£ 


{gd - gu){gd - 


-gi) 


1 -gw{gu + gd] 


+ gugd 


{ge-gu)ige - 


- gd) 



hi + 



lu - gd){gu - ge)' 



{gd - gu){gd - gi) 



6 + 



(104) 



^ ' {gt- gu) {gi - gd) ' 

Here X„ + X^ + = 1 by construction. Note that if these formulae were applied to the 
3HDM-D, they would yield 6^, 6^ and bj, respectively (cf. Eq. [99l in the current model this 
is recovered in the limit ag 0). In particular, the values of all three Xj would necessarily 
lie between zero and one. However, in part of the parameter space of the 3HDM-D+S, one 
of the Xi can be negative. In this part of the parameter space, if one were to (incorrectly) 
assume the 3HDM-D and attempt to solve for the bi, Eq. [99] would fail to yield a solution. 
Thus we see that the footprint of the 3HDM-D+S in the space of h couplings is larger than 
that of the 3HDM-D, and therefore the model with an additional singlet can be distinguished 
from the 3HDM-D in part of its parameter space. 

A negative value for one of the Xi can occur because exactly one of the three denominators 
in Eq. 11041 is negative. This allows us to obtain a lower limit on when one of the Xi is 
negative. We first define 

{gu - gd){gu - ge)Xu if x„ < o, 
= < {gd- gu){gd - ge)Xd if x^ < o, (105) 
{ge - gu){ge - gd)Xe if Xi < 0. 

Note < y < 1 by construction for parameter points where Y is defined. These expressions 
are entirely determined in terms of the h couplings. The lower limit on is then given by 

For completeness, we give here the relations for the parameters bi and aj in terms of the 
h couplings and ^ = 1 — a^: 

- gw{gd + ge) ' - --i^/^ 



bu 
bd 



gdge 

{gu - gd){gu- ge) . 

i-gw{gu + ge) +guge 
{gd - gu){gd - ge) 
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1/2 



^ ^ \ ^ - gw{gu + gd) + 9ugd '\ 
[ {gi-gu){gt-gd) \ 

an = Kgu, ad = bdCjd, ae = bege- (106) 

Because a.^ cannot be uniquely determined in this model, the model is underconstrained and 
the parameters bi and Oi cannot be uniquely extracted. 

These results can easily be extended to models containing two or more singlets by making 
the replacement 

singlets 

We see that it is not possible to tell whether only one singlet or more than one singlet has 
been added to the 3HDM-D on the basis of h couplings alone. 



C. 3HDM-D plus additional doublet (s) 

We finally consider the consequences of adding an additional Higgs doublet $0 to the 
3HDM-D. The additional doublet carries a vev but does not couple to fermions. The con- 
straint equations become, 

al + al + a',+al = l, bl + b^ + b', + bl = 1. (108) 

The normalized couplings of h to SM particles are given by 

_ _ 0,^ _ Old _ ai . , 

gw = aubu + adbd + aibi + aobo, gu = -r-^ gd=-r-, 9e = -r- (109) 

Ou Od oe 

All four bi parameters and gw can be chosen real and positive, while now g^, (jdi and g^ can 
have any combination of signs; in particular, all three of these couplings can be negative if 
ao&o is big enough to keep gw positive. 

Like the 3HDM-D+S, this model is distinguishable from the 3HDM-D in part of its 
parameter space; the parameters and couplings of the current model reduce to the form 
of the 3HDM-D+S in the limit 60 ~^ 0. However, the footprint of the current model in h 
coupling space is larger than that of the 3HDM-D+S, so that in part of the parameter space 
the presence of the extra doublet can be detected, as we now show. 

We again define Xd and in terms of the h couplings as in Eq. 11041 In terms of 
the underlying model parameters, these can be expressed as 

^ ^ ^2 ^ al + blMe - apboigd + ge) 

{gu-gd){gu-gd 
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V l2 «o + bjguge - aoboigu + ge) 



(gd- gu)igd- ge) 

V u2 , (^l + blgugi- aobQ{gu + gd) .,,„^ 
A£ - 0^ H — — — — . (^iiUj 

[ge - gu)[ge - gd) 

Again, + Xd + = 1 by construction. In the limit 60 ~^ 0, these expressions reduce to 
those for the 3HDM-D+S given in Eq. I104( in that hmit the numerator of the second term 
is just Oq, which must he between zero and one. When 60 7^ 0, however, the numerator of 
the second term can be less than zero or greater than one. 

In the part of parameter space with one negative Xi we again construct the quantity Y 
as given in Eq. 11051 In the 3HDM-D+S, Y provided a lower bound for Oq; in particular 
< y < 1 always. In the current model, however, one can also obtain y < (when Xi 
is negative due to the numerator of the second term being negative) or F > 1 (when X^ is 
negative due to the denominator of the second term being negative and the numerator of 
the second term is sufficiently greater than one). Neither of these possibilities can occur in 
the 3HDM-D+S and they therefore allow the current model to be distinguished in part of 
its parameter space. 

The analysis can easily be extended to the 3HDM-D plus two or more doublets. We have 
seen in the case of the 2HDM-I plus an additional doublet (Sec. IIIIPp and the 2HDM-II 
plus additional doublets (Sec. IIV C|) that, once the model already contains one doublet that 
does not couple to fermions, adding additional doublets that do not couple to fermions does 
not change the model footprint in h coupling space. The same is true for the 3HDM-D 
plus additional doublets. Adding one additional doublet changes the model footprint as we 
have seen. Adding a second additional doublet, however, does not further change the model 
footprint; thus it is not possible to tell how many additional doublets have been added to 
the 3HDM-D based only on the h couplings. 

The addition of singlet (s) to the 3HDM-D plus a doublet can be dealt with in a similar 
way. As far as the couplings of h are concerned, adding a singlet is indistinguishable from 
adding an additional Higgs doublet with zero vev. We see then that it is not possible to 
tell whether singlets have been added to the 3HDM-D plus a doublet based only on the h 
couplings. 
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VI. RADIATIVE CORRECTIONS 



In order to translate between the tree-level Lagrangian parameters gw, Qu-, Qd-, and ge, and 
experimentally observable h production cross sections and decay partial widths, radiative 
corrections must be included. This program has been carried out in great detail for the 
SM Higgs as well as for the MSSM. For more general multi-Higgs-doublet models, however, 
detailed results are lacking; such corrections would be needed for a translation between 
observables and the underlying Lagrangian parameters at the few-percent level. We can 
however make the following general observations. 

QCD corrections are universal and can be taken over from the SM, assuming that no new 
strongly-interacting particles contribute. In the MSSM, for example, squarks and gluinos 
yield large flavor-specific radiative corrections; integrating out these contributions into an 
effective Lagrangian yields the formalism but leads to a violation of the underlying 
natural flavor conservation of the MSSM Higgs sector. 

Electroweak radiative corrections are not universal and in principle must be computed for 
each model. These depend on the model content - both the extended Higgs sector and any 
additional new physics that may be present. Some parts of these corrections can be simply 
absorbed into our parameterization; for example, the largest electroweak corrections to the 
MSSM Higgs sector from top quark and top squark loops can be absorbed into an effective 
Higgs sector mixing angle a. However, vertex corrections remain an issue for precision 
parameter extraction. 

Experimental determination of the relative signs of h couplings is also potentially prob- 
lematic. These signs are accessible only through the interference of competing amplitudes 
in loops, and thus nonstandard sign combinations can be masked or faked by additional new 
contributions to loop-induced couplings. 

Throughout we choose the phase of h such that gw is positive. The sign of gu is then 
accessible through the /i77 coupling: in the SM, the W loop dominates while the top quark 
loop interferes destructively, reducing the /i — > 77 partial width by ~ 30%. The relative sign 
of gd and g^ is in principle accessible from the ggh coupling: again in the SM the top quark 
loop dominates, while top-bottom interference is about a 10% effect for moderate Higgs 
masses; however, the QCD scale uncertainty is still of this order. The sign of gp^ will be 
even more difficult, since its contribution to h^^ is extremely small. The loop-induced h^Z 
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coupling would provide additional information, but its experimental detection does not seem 
feasible at the moment. The best strategy may be to examine all possible sign combinations 
for h couplings and enumerate their implications for the underlying model parameters and 
the size of the possible new physics contributions to loop-induced h couplings. 

VII. DISCUSSION AND CONCLUSIONS 

Our ultimate aim in this work is to provide a framework for distinguishing among com- 
peting models for the Higgs sector. To this end we have studied the patterns of tree-level 
couplings of a single CP-even state h in all models that can be constructed out of SU(2)i: 
doublets and/or singlets, subject to the requirement of natural flavor conservation. Distin- 
guishing one model from another relies not only on the underlying theoretical distinctions 
between the values taken by the h couplings, but also on the experimental and theoreti- 
cal precision with which those couplings can be measured. Here we collect our theoretical 
results, then turn to the question of model discrimination based on experimental data. 

Our theoretical results can be conveniently summarized in the form of a decision tree, 
as follows. We assume a deviation from the SM; n > 1 counts additional singlets (S) and 
m > 1 counts additional doublets (D) that do not couple to fermions. We denote the pattern 
relation involving fermion couplings gi and cjj as Pij = gwiOi + 9j) ~ Qidj- The Xi factors 
are deflned in Eq. 11041 and Y is deflned in Eq. 11051 

(i) 9u = 9d = 9e (Type-I-like) 

(a) gw = gf. SM+nS; 2HDM-I when ($o) = 0. 

(b) gw ^ gf 2HDM-I; 2HDM-I+nS; 2HDM-I+mD; 2HDM-I+nS+mD. 

(ii) 9d = 9i^ 9u (Type-II-like) 

(a) P„rf = Pui = l: 2HDM-II. 

(b) 0<Pud = Put < 1: 2HDM-II+nS; 2HDM-II+mD; 2HDM-II+nS+mD. 

(c) P„rf = > 1 or Pud = Pui < 0: 2HDM-II+mD; 2HDM-II+mD+nS. 

(iii) gu = ge7^ gd (flipped 2HDM-like) 

(a) Pud = Pm = 1: flipped 2HDM. 
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(b) < Pud = Pm < 1: flipped 2HDM+nS; flipped 2HDM+mD; flipped 
2HDM+nS+mD. 

(c) Pud = Pm> 1 or Pud = Pm < 0: flipped 2HDM+mD; flipped 2HDM+mD+nS. 

(iv) (ju = Od 7^ 91 (lepton-speciflc 2HDM-like) 

(a) Pui = Pdi = 1: lepton-speciflc 2HDM. 

(b) < Put = Pdi < 1: lepton-speciflc 2HDM+nS; lepton-speciflc 2HDM+mD; 
lepton-speciflc 2HDM-|-nS-|-mD. 

(c) Pui = Pd£ > 1 or Pui = Pdi < 0: lepton-speciflc 2HDM-|-mD; lepton-speciflc 
2HDM+mD+nS. 

(v) gu^9d^ ge 

(a) Pui = l: MSSM with A,,. 

(b) Pui ^ 1 

i. < X, < 1: 3HDM-D; SHDM-D+nS; 3HDM-D+mD; 3HDM-D+nS+mD. 

ii. One of < and < F < 1: 3HDM-D+raS; 3HDM-D+mD; 3HDM- 
D-hnS-hmD. 

iii. One of X, < and r < or F > 1: 3HDM-D+mD; 3HDM-D+mD+nS. 

In particular, we count 15 models (or sets of models) that are distinguishable in principle 
based on the couplings of h. Explicit formulae for h partial widths (equivalently couplings 
squared) in the decoupling limit, i.e., for small deviations from the SM predictions, are 
collected in Table [H 

We now give a flrst comparison of our theoretical results to anticipated LHC and ILC 
measurements of Higgs couplings-squared. In Fig. H] we illustrate the behavior of the partial 
widths (equivalently couplings squared) of /i to or Z boson pairs, down-type quarks, 
up-type quarks, and charged leptons, normalized to their SM values, as a function of the 
decoupling parameter b. Because we consider the full range — 1 < 5 < 1, we use the exact 
formulae from the text rather than the decoupling limit approximations of Tabled We show 
results for the SM plus a singlet (Eq. [25D, the Type-I 2HDM (Eqs. El and [35]), the Type-II 
2HDM (Eq. EHD, the flipped and lepton-speciflc 2HDMs, and the democratic 3HDM. In all 
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Model 




^ d/-^ d 


ul u 




SM 


1 


1 


1 


1 


SM+S 


1 - (5^ 


1 - 52 


1 - (5^ 


1 - 52 


2HDM-I 


1-6^ 


1 + 26/tf3 


1 + 25 /tp 


1 + 25 /tp 


2HDM-II 


1 - (5^ 


1 - 2tp6 


1 + 25 /tp 


1 - 2tp5 


2HDM-II+S 


1 - (5^ - 


1 - 2t^5 - 


1 + 25/t^ - e2 


1 - 2t^5 - 


2HDM-II+D 


1 - (5^ 1 


- 25{s^tf3/cn + 


c^tn) 1 + 25{sy/cntf3 - Cytn) 1 


- 25{s^tp/cn + c^tc) 


Flipped 2HDM 


1 t-9 

1 - (5^ 


1 - 21/36 


1 + 25/t(3 


1 + 2(5/t^ 


Lepton-specinc 2HDM 


1 - (5^ 


1 + 2(5/t/3 


1 + 2(5/t/3 


1 - 2tfj5 


MSSM 


1-52 


1 - 2t'ij6 


1 + 25/t(3 


1 - 2*0(5 


3HDM-D 


1 - ^ 


- 26{s^tf3/cn + 


c^tn) 1 + 25{sj/cntf3 - c^tn) 


1 + 25c^/tQ 



TABLE I: Behavior of the Higgs partial widths (equivalently couplings squared) near the decoupling 
limit, \5\ <C 1. For the 2HDM-II+S we also require <C 1. The other parameters are defined as 
tp = tan /5 = Vf/vQ in the 2HDM-I, Vu/vd in the 2HDM-II, flipped 2HDM, and 3HDM-D, and Vq/vi 
in the lepton-specific 2HDM. For the MSSM we define t'^ = tan/3' = v^il - cot'^ (3 Ab)/vd{l + Af,). 
For the 2HDM-II+D and 3HDM-D we also define cq = cosfl = \Jvl + vj/vsM and 7 is the 
remaining mixing angle that parameterizes the state h. 



models except the SM+S we set tan (3 = 5; for the democratic 3HDM we also set sini7 = 0.2 
(corresponding to Vi = 50 GeV) and cos 7 = 0.5. 

On the right-hand side of each plot in Fig. H] we also show the expected la measurement 



uncertainties of the squared Higgs couplings at the LHC from Refs. ll|, |87|] (summarized in 
Table [Xljl , assuming SM coupling strengths and a Higgs mass of 120 GeV. The coupling fit 
of Refs. P, was based on anticipated Higgs production and decay rate measurements 
from the LHC using 300 fb"^ of integrated luminosity times two detectors. Vector boson 
fusion channels (which have only been studied for 30 fb~^ to date) are scaled to 100 fb~^ 
to account for potential degradation at high luminosity running. The fit assumes SM rates 
in all channels and allows for an unobserved component of the Higgs total width as well as 
nonstandard contributions to the ggh and /i77 vertices. It further assumes = (jz — 1-05, 
which is valid for models containing only doublets and/or singlets. Theoretical uncertainties 
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FIG. 4: Higgs partial widths (equivalently couplings squared) as a function of the parameter b in the 

SM+singlet, the Type-I, Type-II, flipped and lepton-specific 2HDMs, and the democratic 3HDM. 

We set tan/3 = 5 for all models except the SM+S; for the 3HDM-D we also set sinfi = 0.2 and 

cos 7 = 0.5. At the right of each panel we show the expected Icr LHC measurement uncertainties 

for the Higgs couplings-squared to WW (blue), hh (red), ti (green), and rr (orange) for ra^ = 120 

GeV and SM event rates, taken from Table HH Note the log scale on the y axis. 
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9w al 



9t 



9r 



LHC [11] 22% 43% 32% 27% 
ILC [161 2.4% 4.4% 6.0% 6.6% 



TABLE II: Expected uncertainties on Higgs coupling-squared measurements at the LHC and ILC, 
assuming = 120 GeV and SM rates for all processes involved. See text for details. 

on Higgs production rates due to QCD scale uncertainty were also included. 



The results of Refs. 



11 



87| will likely change when updated experimental and theoretical 



results are included. Updates of all experimental channels are now available in the CMS 



Physics TDR 



and the ATLAS Computing System Commissioning (CSC) Notes [89i]- In 



particular, the critical tth, h bb channel is dead; more work is needed on the experimental 



side to evaluate the potential of newly-proposed bb final state channels like Wh, h bb 90|. 
Progress has also been made on the higher-order corrections to the gg —>■ h cross section 91 1. 
A new fit involving these more sophisticated results will require some work. 

Nevertheless, we sketch the current situation as follows. We scan over model parameters 
and compute a relative to the SM prediction according to 

i=W,b,t.T 



X 



[ST 



;iii) 



using the LHC uncertainties in the partial widths from Refs. [111 . I87l | as summarized in 
Table Hn We make no attempt here to account for the correlations in the extracted couplings. 
One-, two- and three-sigma contours are shown in Fig. [5] for the LHC. For comparison, we 
show the corresponding ILC expectations in Fig. [6l 

The SM plus a singlet contains only one additional parameter that universally shifts the 
partial widths to all SM decay modes, while the 2HDM models listed contain 6 and tan (3 as 
free parameters that describe the Higgs coupling. Since the 3HDM-D model has four free 
parameters, 6, tan/5, Q, and 7, we marginalize over Q and 7 by evaluating a Markov Chain 
Monte Carlo (MCMC) [921, MM following the procedure of Ref. |94|. 

In summary, we have provided a first roadmap to determine the underlying model of 
electroweak symmetry breaking under the assumption that only Higgs doublets and/or sin- 
glets participate and the Glashow-Weinberg-Paschos condition for natural flavor conserva- 
tion holds. Our approach is based on the couplings of a single identified CP-even Higgs 
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5 10 15 20 5 10 15 20 

tanp tanp 



FIG. 5: Regions of parameter space with combined h couplings within one, two and three a 
corresponding to the inner, middle and outer contours, respectively, of the SM limit for various 
models, based on the expected LHC sensitivities given in Table |TI1 Values of are calculated 
according to Eg. Illli The 3HDM-D model contains four free parameters, 5, tan/3, fi, and 7; we 
marginalize over Q and 7 by evaluating a Markov Chain Monte Carlo (MCMC). The minor wiggles 
in the shapes of the contours is due to the numerical precision of the MCMC. 
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FIG. 6: Same as Fig. [5]but for the ILC, using the precisions on coupHngs squared given in Table Ull 



state without regard to other Higgs particles that may appear in the spectrum. We restrict 
our considerations to tree-level decays of the Higgs boson to avoid complications from new 
physics that may appear in loop-mediated decays. We described 15 classes of models and 
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compared their predictions for the shifts in the Higgs couphngs relative to the SM. In each 
case, we presented formulae for the couplings of a single CP-even state htoWoiZ boson 
pairs, up-type quarks, down-type quarks, and charged leptons as a function of the model 
parameters at tree level. Where possible, we also inverted those relations to provide explicit 
formulae for the model parameters in terms of the h couplings. We summarized our results 
in a decision tree that can be used to differentiate among the models. While extraction 
of the couplings of h with sufficient precision at the LHC will be challenging, our results 
provide a starting point for a more detailed study of model discrimination based on future 
experimental results. 
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